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Abstract 

A numerical study of coupled to the dilaton field, static, spherically symmetric 
monopole solutions inspired by the Kaluza-Klein theory with large extra dimen- 
sions are presented. The generalized Prasad-Sommerfield solution is obtained. We 
show that monopole may have also the dilaton cloud configurations. 



1 Introduction 

Recently there has been considerable interest in the field theories with large extra 
spacetime dimensions. In comparison to the standard Kaluza-Klein theory these 
extra dimensions may be restricted only to the gravity sector of theory while the 
Standard Model (SM) fields are assumed to be localized on the 4-dimensional 
spacetime (Antoniadis, Arkani-Hamed, Dimopoulos and Dvali (1998); Antoniadis, 
Arkani-Hamed, Dimopoulos and Dvali (1998); (Arkani-Hamed Dimopoulos and 
Dvali (1999)). It is promising scenario from the phenomenological point of view 
because it shifts the energy scale of unification from 10 19 GeV o 10 — 100 TeV. 

The gauge field theory is extended by inclusion of the dilatonic field in such 
of theories. Such fields appear also in a natural way in Kaluza-Klein theories 
(Appelquist, Chodos and Freund, (1987)), superstring inspired theories (Witten 
(1985); Ferrata, Lust and Teisen (1989)) and in theories based on the noncommu- 
tative geometry approach (Chamseddine and Frohlich (1993) ). 

As previous studies have already shown the inclusion of a dilaton in a pure 
Yang - Mills theory has consequences already at the classical level. In particular 
the dilaton Yang - Mills theories possess 'particle - like' solutions with finite energy 



1 



which are absent in the pure Yang - Mills case. Analogous equations have recently 
been obtained for the 't Hooft - Polyakov monopole model coupled to the dilatonic 
field (Lavreashvili and Maison (1992)). 



2 The dilatonic gauge field theory 

Dilatons appear in the higher dimensional theory after the process of the spacetime 
compactification. The main idea of the theory with large extra dimensions is that 
the gravity is realized in the more dimensional spacetime (the bulk) while the 
matter is confined in the four-dimemsional spacetime (the brane). To be clear and 
simple, we will consider the six-dimensional gravity. Let us now consider the action 
integral of Einstein- Yang-Mills-Higgs theory in the six-dimensional spacetime: 



J d 6 xV=g~ 6 L, (l) 



where ga = det(gMN) and M = {/x, i}, N = {u, j} with x M = i = 1,2. 

The metrical tensor in the six-dimensional spacetime can be written: 

9MN ~ { _ r 2^. e +2^)//o ) > ( 2 ) 

According to above definition we can write: 

v^=^r 2 2 e- 2 ^. (3) 

In the equation (||) 

9,u = e- 2fW//o V ( 4 ) 
represents the four-dimensional metric in the Jordan frame while ~g^ v in the Ein- 
stein frame. We consider here the Lagrangian of the Einstein- Yang-Mills-Higgs 
field as follows: 

L = L g + L Y MH${y), (5) 

^ = -i(«- 2A )> (6) 
Lymh = -\F« v F a » v + \D^ a D^ a - £/($), (7) 

where kq is the six-dimensional gravitational coupling. L, L g , Lymh - describe the 
total Lagrange function, the gravity in six-dimensional spacetime and the Yang- 
Mills-Higgs field parts on the brane emdeded in six-dimensional space, respectively. 
In general non-vanishing cosmological constant (A ^ 0) is possible. This case leads 
to the interesting monopole solution (Lugo and Shaposhnik (1999, 2000)). In our 
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paper we shall focus our attention on the A = case. All calculations should 
include metric guN ©, so for example = g^ v D v . Let us compactify the 
six-dimensional spacetime to the four-dimensional Minkowski one on the torus 
(Mq — > M4, x S l x S 1 ). In this paper we assume that the extra dimensions are 
compactified to two-dimensional torus with a single radius r 2 . The six-dimensional 
action may be rewritten as: 

S = Jd 4 xJ d 2 y^g- 6 L = J d 4 x^C, (8) 

where J d 2 y = (27rr 2 ) 2 and C is the effective Lagrange function in four-dimensional 
spacetime. The six-dimensional gravitational coupling kq = 8ttGq is convenient to 
define as 

where M - is the energy scale of the compactification (~ 10 — 100 TeV). Com- 
pactiflcation of the six-dimensional gravity on the torus gives the Lagrangian (||) 
for the four-dimensional gravity 

in the Einstein frame, where 

1 (2vrr 2 ) 2 



(10) 

K HQ 

is the four-dimensional coupling constant or k = 8ttG n = 8vrM~f . From the 
above equation ([R]) we get 

M 2 Pl = A-nM 4 rl. (11) 

Cosmological consideration (Hall and Smith (1999)) gives the bound M ~ 100 TeV 
which corresponds r 2 ~ 5.1 x 10~ 5 mm from equation (11). Compactification 



of gravity on the five-dimensional spacetime is rather unpysical (r 2 ~ 10£;m), 
however the nice five-dimensional spacetime compactivitation was propsed recently 
(Randal and Sundram ( 1999a, b)). The Planck mass Mpi ( |Tl|) is not longer a 
fundamental constant, it may change during the Universe evolution (Flanagen Tye 
and Wasserman (1999)). For four-dimensional Minkowski spacetime {jg^ v = Tj^u) 



%) = ^ {x)/f °{-d^i + (12) 



The last term in the equation (12) can by transform into the first one by differen- 
tiating by parts. The parameter /o (or re-scaling of the £(x) field) is determined 
by the Planck mass (at present time) as: 
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fo = ^=M Pl ~ 4.87 10 18 GeV/c 2 (13) 



to produce the 1/2 term in for the dilaton field in (16). The /o parameter deter- 
mines the dilaton scale fo- At the present time fo is rather high, so the interaction 
with dilatons can be neglected. However, in the early universe when the Planck 
mass Mpi was smaller (for details see (Flanagen Tye and Wasserman (1999))) also 
the value of the fo was smaller, too. 

As a result of compactification of the six-dimensional Lagrangian we get the 
Lagrange function for the Yang-Mills-Higgs fields. Fluctuations around the four- 
dimensional Minkowski 

9^=V^ + h ^u(x,y) (14) 
will produce interaction with Kaluza-Klein dilatons 

n 

with the typical mass scale M (for n l ^ 0). 

In this paper we shall apply this approach to the simplest SO (3) gauge field 
theory. The 50(3) gauge field theory has nice monopole solutions ('t Hooft (1976)) 
which have produced trouble in the cosmology and have been the reason to intro- 
duce the idea of inflation. The main idea of this paper is to examine how the 
monopole solution looks like in the dilatonic gauge field theory inspired by the 
Kaluza-Klein gravity with the TeV scale. 

The dilatonic gauge field theory may described by the Lagrangian (defined 
by equation (§), in the first approximation when g^ = rj^ in the Minkowski 
spacetime) we get 

C = \d^{x)d^{x) - y 2 ^)/foF^F a ^+ (16) 

where for 50(3) theory we have: 

U(&) = ^($ a <S> a -v 2 ) 2 (17) 

with the 50(3) field strength tensor F^ u = d^W? - d v W% + ge abc W b W^. The 
50(3) gauge symmetry rotates the Higgs field $ a . The covariant derivative is 
given by D^ a = d^ a - ge a b c W^ c . Now we have = rf v D v . The Higgs 
potential has degenerate true vacua forming the sphere S 2 ($ a $ a = v 2 ). The 
Euler-Lagrange equations for the Lagrangian (9) are scale - invariant: 
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x v _> x >» = e ^x», (18) 

Z^? = Z + u, (19) 
$ $' = $, (20) 

W« - W£ = e"£w£. (21) 
These transformations change the Lagrange function in the following way: 

_'2u 

C -► £' = e /o£. (22) 

This symmetry can be formulated equivalently as a scaling symmetry on the co- 
ordinates, and the dilaton is often denoted as a Goldstone boson for dilatations. 
The origin of the symmetry of the equations of motion is easy understood from 
the Kaluza-Klein origin of the action. The scale transformations are equivalent to 
a rescaling of the internal dimensions. 

3 The dilatonic monopole. 

The monopole scalar field configuration: 



= v h(r)n a = vH{r) — , (23) 
gr 



where 



n a = — (24) 
r 

describes the 'hedgehog' structure n a and scalar spherically symmetric field H(r). 
The 50(3) gauge field is described by the K(r) field: 

A1 = -e aij —n\l-K{r)). (25) 
gr 

The dilaton field is described by the S(x) function: 

£{x) = f S(x) (26) 
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(if we introduce the dimensionless variable x = gvr). The field Euler-Lagrange 



equations generated by Lagrange function (16) for H(x) and K(x) and S(x) are: 



H"(x) - ^H(x)K(x) 2 + 4e~ 25(r) (^ 2 - H 2 (x))H(x) = 0, 



(27) 



K"{x) + 2K'{x)S'{x) - \e- 2Six) H(x) 2 K(x) - ^^(K(x) 2 - 1) = 0, (28) 



S"(x) + -S'{x)-a- 1 \e 2S{ - x) {{l-K{x) 2 ) 2 + 2x 2 K\x) 2 } (29) 



+ —sol e 



l„-2S(r) 



(1 



H 2 (x), 2 



0. 



In the dilatonic monopole we have two independent dimensionless constants: 

_ A_ 



« = (-) • 

The mass (or the lowest energy) of the monopole in the rest frame is: 

Attv 



Mr, 



with the energy density given by: 



9 



p(x)x dx, 



p{x) = l -aS'(x) 2 + ^e 2S ^{{l-K{x) 2 ) 2 + 2x 2 K'(x) 2 } 
+ ±- A {2H(x) 2 K{x) 2 + (xH'(x) - H(x)) 2 } 
+ * -2S(x) (1 _£!M)2. 



(30) 
(31) 

(32) 



(33) 



Inside the monopole (according to the equations (27), (28) and (pS|)) the asymp- 
totical behaviour when x — > is given as: 



K(x) = 1 -tx 2 + 0(3), 



(34) 



ff(ar) = -ux 2 + 0(3), 



(35) 
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S(x) =a + bx 2 + 0(3), (36) 
where u, t, a are local parameters and b must be determined as: 

2Ae 2a t 2 - ee~ 2a 



12a 



(37) 



Far from the monopole core, if r — > oo (x — > oo), both functions H(x) and K(x) 
should describe the normal vacuum ($ a $ a = v 2 ) with H — ► x and K — ► according 
to the if, X function definition (23|, 25) remembering that x = gur. In this limit 



the energy density rt33n has a simple limit 



P (x)= 1 -aS'(x) 2 + ^e 2S ^. 

The monopole mass in this limit may be rewritten as: 

M mon = lffp(x)x 2 dx = 
*f KV^xS' + le s ) 2 dx + ^^(e 5(0) - e s (°°)). 
The first term vanishes if the dilaton field obeys the Bogomolny equation 

yfaxS' + -e s = 0. 

x 

This equation has a nice solution in the uniform normal vacuum (Bizon (1993)) 

S b = -ln{{^e- s{oo) - l/x)/y/a). (38) 

When r — > oo dilaton field should disappear in the true vacuum. This demand 
gives S(oo) = 0. So, When r — ► oo we have the asymptotic behaviour of the 
solutions 

H(x) = x{l-we- VTex ) + 0{l/x), (39) 
K(x) = ze~ x + 0(l/x), (40) 

S(x) = -ln(^e~ sico) - -^-) + 0(l/x). (41) 

y/ax 

Even when S(oo) ^ it may be removed by the dilaton transformation (|l9|). We 



may solve the differential equations (|27j, 28, |29|) by the iteration method expanding 
them with respect to e: 

oo 

H(x) = Y / z n H n (x), (42) 

n=0 
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K{x) = Y j e n K n {x), (43) 

ra=0 
oo 

S(x) = a- 1 Y J z n Sn{x). (44) 

n=0 

In the first step (n = 0) we obtain the equations: 

H'£{x) - ^H (x)K (x) 2 = 0, (45) 



KS(x) + lK' (x)S' (x) - ^e- 2S ^/ a H (x) 2 K (x) (46) 
-i^l(K (x) 2 -l) = 0, 

9 p 2S (a;)/Q 

S// (s) + -S' (x) ^{(1 " ^o(^) 2 ) 2 + 2x 2 K^(x) 2 } = 0, (47) 

leading to the Prasad-Sommerfield solution (without the dilaton field S(x)). When 
a — > oo we have the Prasad-Sommerfield solution (Prasad and Sommerfield (1975)) 
and we can find for Hq(x) and Kq(x) easily as : 

H (x) = x/tanh(x) - 1, (48) 



Kq(x) = x/sinh(x). (49) 
Finding a nice analytical solution for the dilaton field (Fig. 3, dotted line) 

S (x) = a+S^ + —^^ (-i + 2x 2 + 2xQ D + cosh(2x) (50) 
—2xQd cosh(2x) — 2xsinh(2x)) 

was a crucial point of this paper. The leading term for the dilaton field at infinity 
will be a Coulomb one 

S (x) = — + 0(l/x) 

X 

Qd is the dilatonic charge which originates from the global scale transformation 
(|18-|21). The similarity is striking, but we should remember that an electric charge 
comes from the exp(iea) 6 U(l) gauge symmetry. The global scale transformation 
(18-21) is generated by the exponetrial transformation exp(u//o). However, the 
asymptotic bevaviour at x — > (j3q) admits Qd = 0. 

We present the numerical solutions of the coupled set of differential equations 
( |27| . 28, |29| ) in the next section. 
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4 Numerical solutions 



To solve the monopole equations numerically, we need the starting point (Press, 
Teukolsky and Vetterlino (1992)). To find the starting conditions we can use the 
solutions found from the variational procedures or from the Prasad-Sommerfield 
approximation (|48|, |4^, 50). The trial solutions depending on the variational pa- 
rameters must be postulated in such a way to fulfill the boundary conditions close 
to the center ( p^j35|j36| ) of the monopole and at far outside (|^J4^,41). We postu- 



late the trial solutions: 

H v (x) = x (nx + x\l -e-^)) ^ ^ + 
(1 + x z ) 

K v {x) = (1 ~+lt? 4) ~ 1 " ^ + ° (3) (52) 

and 

(a + bx 2 + Q D x 5 ) , h 2, nf o\ /m\ 

bJx) = 777 ~a + 6x + 0(3), (53) 

(1 + x°) 

where 0(3) are corrections of the third order. For these functions the monopole 
mass was calculated and the trial functions with minimal energy was found. For 
the monopole without dilatons we get the M mon mass of monopole (if S(x) = 0): 

M mon = 15.886 x 10 15 GeV. 

The trial functions give the dilaton configuration close to the monopole case with- 
out dilatons. The minimal variational configurations for such a dilatonic monopole 
for Qd = and v = 10 16 GeV are presented on the Table 1. This shows that the 
mass of the dilatonic monopole dependences on the parameter a (/o) and reaches 
the local minimum (for a ~ 1) lower than for monopole without dilatons. 

For the dilatonic monopole the numerical method was independently verified 
using the Chebyshev polynomial expansion (Michaila (1999)). The monopole so- 
lutions for H (x) and K(x) are known very well so our attention is focused on the 
dilaton solution S{x) especially. The behavior of the H(x) and K(x) solutions 
determined by the boundary conditions is the same as it is in the presence of the 
dilaton field. 

The Chebyshev method allows to calculate the exact solution of the differential 
equations for the discrete set of points. The trial function provides the starting 
data for the numerical solution of the ordinary differential equation (ODE) (the 
shooting method (Press, Teukolsky and Vetterlino (1992) )) or the Chebyshev 
functions method. After this preliminary numerical calculation the method based 
on the Chebyshev polynomial was used. 
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a 


u 


t 


z 


M da (10 15 GeV) 


2.37xl0 5 


0.2396 


0.6021 


2.46903 


16.4125 


2.37xl0 3 


0.2396 


0.6021 


2.46886 


16.4107 


2.37 


0.2591 


0.5769 


2.2637 


15.0284 


1 


0.2879 


0.5389 


1.92221 


14.4402 


0.9 


0.2937 


0.5327 


1.8099 


14.4913 


0.8 


0.3011 


0.52609 


1.7035 


14.6313 



Table 1: The dependence of the monopole mass and the parameters u, t, z 
on the a parameter. 




Figure 1: The dependence of the gauge field K{x) on the x parameter. The 
dotted line corresponds to the Prasad-Sommerfield solution (|49|), the solid 
one to the variational solution flo"2|), the dashed line presents the solution of 
the first iteration. 
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Figure 2: The dependence of the Higgs field H(x)/x on the x parameter. The 
dotted line corresponds to the Prasad-Sommerfield solution (|48|) , the solid 
one to the variational solution (|51|), the dashed line presents the solution of 
the first iteration. 
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Clenshaw and Curtis have proposed almost forty years ago an integration 
method based on the Chebyshev polynomials of the first kind of degree j, 

tj(x) = cos(j * arccos(x)). (54) 

Since then, these methods have become standard. Since the Chebyshev polynomi- 
als are orthogonal and allows to rewrite the function f(x) as 

f( x ) = HVi( x )> ( 55 ) 

3=0 



where (for j = 0) 



and (for j / 0) 



1 " 

ao = - f( x k)to(xk), (56) 
n k=i 



2 n 

a i = -^lf( x k)tj{xk)- (57) 
n k=i 

The grid of n points Xk are zeros of the Chebyshev polynomial tj(x). This decom- 
position allows us to present the derivative of the function f(x) as 

f'(x l ) = Y / D ik f(xk), (58) 
k 



where the matrix 



n— 1 



Dtk = J2 -tjixk^ixi), (59) 



3=0 C 3 

and cq = n, Cj = n/2 (at j / 0). This fact transforms the ordinary differentional 
equation: 

+ p(x)^f- + q(x) = r(x) (60) 
dx £ ax 

into an appropriate linear equation: 

^Aikfixk) = Ui- (61) 
k 

So, we can have exact solution for discrete number of points. This method may 
be also used to nonlinear equation 

d 2 f df 

m d^ +p{x) ic +F{Lx) = Q - (62) 
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If we have a starting function /q then we expand / around /q 



f(x) = f (x)+e(x), (63) 



and approximate the eq.(p2|) with the eq.(|60D and then solve numerically. The 
solution may be treat now as a starting function for the next iteration, and so on. 
The iteration may last so long as an arbitrary precision is reached. 

The perturbation around e produces the series of differential equations 

fa,b + ^2Pab,ni X )fb,n(x) + 1ab,n{x) fb,n( X ) = r a,n{ X )> ( 64 ) 
b b 

where the vector 

/„ = {H n (x), K n (x), S n (x)}. (65) 

For example when n = 1 the first equation (a = 1) corresponds to pib = 0, qu = 
2K 2 (x)/x 2 , q l2 = 2H (x)K (x)/x 2 , r x = exp(-2S (x)/a)(H 2 (x) - x 2 )H (x)/x 2 , 
and so on. 

In the monopole case the starting function are these obtained by the variational 
method. After expanding around trial functions (51,52,53) we obtain system of the 



differential equations of the type (p0|). After that the numerical solution may be 
obtained on the grid of the x^. The numerical solution for the dilaton field found 
by the Chebyshev numerical method is presented on the Fig. 3 (the solid line). 



5 Conclusions 

The aim of this paper was to present a numerical study of the classical monopole 
solutions of the SO (3) theory coupled to the dilaton fields. We have shown that a 
monopole is surrounded by the dilaton cloud S(x). In the field theories with large 
extra spacetime dimensions the Planck mass is not longer a fundamental constant 
and may change itself during the evolution of the universe. As a consequence the 
parameter /q changes, too. We have shown that the dilatonic monopole reaches the 
minimal mass when fo~v with the mass lower a bit than for monopole without 
dilatons. 

There is analytical solution in the Prasad-Sommerfield limit. 

The spherically symmetric dilaton solutions coupled to the gauge field or grav- 
ity are interesting in their own and may moreover influence the monopole catalysis. 

However in the theory inspired by the Kaluza-Klein theory with large extra 
dimensions also new interaction with massive (~ M ) Kaluza-Klein gravitons takes 
place. In the four dimensional spacetime the monopole solutions is stable due to 
the monopole topological charge. Now the interaction with Kaluza-Klein gravitons 
h™{x) may cause disintegration of the monopole. 
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Figure 3: The profile of the dilaton clound S(x). The solid line presents 
the Chebyshev numerical solution S c (x) and the dotted line is the Prasad- 
Sommerfield solution Sq(x) (50). 
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